28/05/08 



. . . The integrable 0(6) model and the correspondence: checks and predictions 

00 
O 

o 

■ Francesco BuccheriQ and Davide Fioravanti0 

OO ■ Sezione INFN di Bologna, Dipartimento di Fisica, Universitd di Bologna, 

^ Via Irnerio 46, Bologna, Italy 

Abstract 

We exactly compute the energy density of the integrable 0(n) non-linear sigma 
, ■ model as a convergent series. This series is specifically analysed for the very 

important 0(6) symmetry, since it was suggested to result as a peculiar limit of 
^ , the AdS string theory by Alday and Maldacena [1]. In this respect, the 0(6) 

! model gives also refined confirmations and predictions once compared with the 

^ ; SYM Bethe Ansatz [H [3]. 

o 

00 

O 

> 

X 



* Laurea student of the University of Bologna " Alma Mater Studiorum" . 
tE-mail: bucca.virus@tin.it and fioravanti@bo.infn.it . 



1 



1 The questio 



The planar s/(2) sector of A/" = 4 SYM contains local composite operators of the form 

Tr(P'Z^) + ...., (1.1) 

where V is the (symmetrised, traceless) covariant derivative acting in all possible ways 
on the L bosonic fields Z. The spin of these operators is s and L is the so-called 'twist'. 
Moreover, this sector would be described - via the AdS/CFT correspondence |4] - by 
rotating string states on the AdSs x spacetime with AdSs and charges s and L, 
respectively [5l[6]. Proper superpositions of the operators (11. ip have definite conformal 
dimension A depending on 't Hooft coupling A = Stt'^q'^ 

A = L + s + ^{g,s,L), (1.2) 

with 'j{g, s, L) the anomalous part. In fact, the correspondence would assign to this 
dimension A the exact value of the energy density of a rotating string, provided the 't 
Hooft coupling is identified with the string tension: -\/A = This is indeed a duality 
relation between the coupling constants involving the semi-classical expansion on the 
string side O [6] . 

A great boost in the evaluation of the anomalous 'y{g, s, L) has come from the 
discovery of integrability and thus of a bethe Ansatz, although in another sector, the 
purely bosonic so(6), and at one loop of the gauge theory [7j. On the other hand, in the 
twist sector of one loop QCD the (integrable) Bethe Ansatz problem was at hand [SI E] 
and later on it has been realised to be equivalent to its supersymmetric relative with 
the occurence of the integrability extended to the whole theory and at all loops [TO] . 
This is the scenario on the side of the SYM theory in the sense that, for instance, any 
operator of the form (11. 11) is associated to one solution of some {asymptotic^ Bethe 
Ansatz-like equations and then any anomalous dimension is expressed in terms of this 
solution. 

As a confirmation of the correspondence, integrability has been also uncovered and 
studied in the superstring theory [12] and in this respect our interest will be limited 
to the semi-classical calculations. In this approach, the string tension diverges since it 
plays the role of the inverse of the quantum Planck constant. Therefore, the A +oo 
limit yields a power expansion in 1/ a/A [13] . Which, in particular means, that it needs to 
be implemented before any other limit and thus endowing the semiclassical calculations 
with a different limit order with respect to the gauge theory (cf. below for more details 
and [21 [3|). 

On both the string and the gauge theory, an important double scaling may be 
considered: ^ 

s ^ oo , L ^ oo , j = — — = fixed . (1.3) 

2 In s 

^This important limitation emerged because the Bethe Ansatz is reahsed by somehow using the 
on-shell S-matrix [llj . 
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In fact, the relevance of this logarithmic scaling for the anomalous dimension has been 
pointed out and deeply studied in [13] and [1] within the semi-classical string theory 
(cf. also [Hj within the one-loop SYM theory). Moreover, these long operators with 
large spin have been recently shown to satisfy the Sudakov scaling [H [15] 

7(^?, s, L) = f{g, j) In s + 0((ln s)"^) , (1.4) 

which generalises the one loop result of [H]. Actually, in [15] this statement was ar- 
gued by computing iteratively the solution of some integral equations and then, thereof, 
the generalised scaling function, f{g,j) at the first orders in j and g^: more precisely the 
first orders in g"^ have been computed for the first generalised scaling functions fn{g), 
forming 

oo 
n=0 

As a by-product, the reasonable conjecture has been put forward that the two-variable 
function f{g,j) should be analytic (in g for fixed j and in j for fixed g). In [18] similar 
results have been derived for what concerns the contribution beyond the leading scaling 
function f{g) = fo{g), but with a modification which has allowed to neglect the non- 
linearity for finite L and to end-up with one linear integral equation (LIE). The latter 
does not differ from the BES one (which covers the case j = 0, cf. the last of [ID]), but 
for the inhomogeneous term. Moreover, a suitable modification of this LIE has been 
applied in [2l[3] to derive still a LIE in the scaling (11. 3p (for any g and j). This is indeed, 
a way to determine the generalised scaling function f{g,j) and also its constituents fn{g) 
for all values of j and g, thus interpolating from weak to strong coupling. Today, an 
interesting paper [16] appears which seems to have some equivalent equation, coming 
from [15], from which it apparently derives a map to the 0(6) sigma model and the 
leading strong coupling behaviour of /s (as /2 = appeared already in [15] and fi in 

my 

In the following, we will constrain ourselves to the analysis of the 0{N) energy 
density of the non linear sigma models, since one representative, the 0(6) model, was 
suggested by Alday and Maldacena [1] to represent the limit theory with small 5*0(6) 
charge j ^ V^- For in this regime the masses of the fermions do not contribute to the 
energy density, but to fo{g), where they give the natural UV cut-off. Therefore, we will 
show that the additional energy density, 2Q{g,j) (cf. (18.751) ) contained in f{g,j) can be 
computed exactly at least at the leading mass gap [m of 12. 71) order by means of the 0(6) 
computations. Hence, we will compute below the 0(6) energy density as a convergent 
series in j/m, checking a perfect agreement with the gauge theory computations of [3] 
up to the first interaction and model depending term f^ig), i.e. An all order explicit 
match would be desirable for the future and in this respect the following 0(6) model 
results give a series of exact predictions. 

^0{\ns~°°) means a remainder which goes faster that any power of Ins: Um (In s)'^0((ln s)~°°) = 

s — >oo 

0,Vfc > 0. 
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2 The 0(N) nonlinear sigma model 



We want to show that the result in [TH] for the energy density f2 as a function of the 
particle density j in the strong coupling regime {g » 1) with j << m, namely the 
"non-interacting fermion gas" approximation, keeps correct when all orders in B are 
retained in the calculation, i.e. that it is the correct result up to the third order j^. In 
formulas, that 

where the mass gap is (see [21] and [1], with = j) 

93/47^1/4 , 

We begin in full generality, by considering a 0(N) nonlinear sigma model. A manip- 
ulation of BA equations [19] provides a " pseudoenergy" e{'d) as the solution of a linear 
integral equation of Fredholm type 

ei-d) - [ Ki-d - ^')e{'&')d^' = h - m cosh 

J-B 

With the boundary condition 

e{B) = (2.8) 

The condition (12.81) allows us, at least in principle, to determine the value of the pa- 
rameter B. The kernel K {i}) comes from the Eg two-particle S-matriJ^ 

..0, r(i + f)r(i-f)r(i + A + f)r(A-f) 

^ ^ r(i-f)r(i + f)r(i + A-f)r(A + f) ^ ''^ 

through 

K{^) = ^^logSi^) (2.10) 

where A = What is more, the Bethe Ansatz procedure gives the Helmholtz free 
energy / as 



— f(h,m)= [ cosh^e(^)d^ (2.11) 
m J 



-B 



while general Thermodynamics provides the density of particles 



J = -^f{h,m) (2.12) 



^This matrix contains an ambiguity of the ODD kind [17], as underlined by Zamolodchikov and 
Zamolodchikov. What people generally do, is to use the most simple solution (j2.9p 
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Gathering up all the elments, an optimistic procedure would be that of solving ([2]) for 
e, put it into (12. lip , then calculate j by (12.121) and, finally, Q by Legendre transform 



n{j,m)=f{h,m)+jh (2.13) 

The "non-interacting fermion gas" approximation corresponds to the limit 5 — > in 
(E]), i.e. ^ ^ 0+. Indeed, considering S{^) = -1 and K(i9) = s^'^^i'^) = 0' 
equvalent to saying that B ^ allows one to retain only the forcing therm of ([21), that 
is the zero-th order approximation of the Liouville-Neumann procedure. 



3 The ID non-interacting fermion gas 

Despite its simplicity, this model is noteworthy because it represents a paradigmatic 
example of how Bethe Ansatz works and also because it produces results that can di- 
rectly be compared with those coming from more advanced calculations from the SYM4 
front. What we will call "fermions" are introduced in a simple fashion by defyining 
their two-particle S-matrix as 

S{^) = -1 (3.14) 

now we can set in motion the calculation sketched above, with K^'d) = for (12.101) . 
Then, by fl2l). fl212ll . (EH]), fl213D we have 



e{-&) = h — m cosh 



(3.15) 



f{h)—f{0) = / (h— m cosh {}) cosh (I'd = [2h sinh B — rriB — m cosh B sinh B] 

2n J 2tt 

(3.16) 

B 

j = — I cosh ■dd'd = — sinh B B = arcsin h ( — j ) (3-17) 



27r 



-B 



Vm" 



n 



2 2 

Tfi TTl 

f + jh = — [B + sinh B cosh B] = — 

2tt 27r 



m 

2^ 

2 

m 
vr 



n l/7r\3 3/7r 

— J - 7 — J ) + TTT — i 
m b \m / 40 Vm 



/ TT 



TT 



arcsin /i — j H J \/ 1 + — J 



vr . / 1 / TT . 

m \ 2 Vm 



1 / TT . 
8 \rn 



vr 1 
m D Vm 



TT \ 3 

3 



1 / TT .\ 5 

40 \m 



O 



m' 



TT 



0(4 

m' 



(3.18) 



The series contains only the odd powers of j: we will see that turning on an interaction 
will turn on the coefficients of the even, powers from the fourth power on, together with 
affecting the odd terms. 
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4 The relativistic interacting gas: some preliminar- 
ies 

The formal solution of ([2]) is easily seen to be 

£(^) = h-m cosh ^ + J K{d - d') [h - m cosh d' + j K{§' - d")e{d")dd''^ d§' = 
= . . . = h — m cosh i) + 

B B 

+ J2 ^^i-- / d^nK {^-^i)K{^i- ^2) ■■ -K i^n-i-^n) {h-m cosh ^n) 

00 \ \ 

= 1 + ^B"" / dxi... / dxnK{Bx - Bxi)K{Bxi - BX2) 



. . . K{Bxn-i - Bxn) {h-m cosh(5x„j()1.20) 



where 
K{^) 

with 



47r2 



(4.21; 



ip {x) 



r (x) d 



r fx) dx 



IriT {x) = -7 + V (— 

^ \n + l n + x J 

n=0 ' 



where 7 is the Euler-Mascheroni constant. Its /;,- derivative, by the aid of ([2]) and fl2.8p 
reads 



dh 



l + ^K{^-B)e{B)) — + ddK{d-d', 



d^K{^-^')^^ 
oh 

-B 



(4.22) 



that is, another Fredholm integral equation, with a constant forcing therm. Its solution 
is 

B B 



dh 



l + Y. I dA... I d^nK{^-^i)K{A-^2)...K{^n-l-^n) 
n=l „ n 



-B ~B 

1 1 



1 + / dxx... \ dxnK{Bx - Bxi)K{Bxi - BX2) 

n=i 

...K{BXn^l-BXn) (4.23) 
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By the use of the formula 



-X 



'iP ( -X 



fc=0 



k + X 



we can rewrite the Kernel as 
K{x) -- 



^ OO 



2tx l^n{k + 1) + - ^Tr{k + 2A) + - 

k — fc — 



+ 



(-1)' 



_'Tr{k + l)-i{^-^') ^Ti{k + 2A)-i{'&-^') 

k — k — 



1 

2^ 



'iy2iT{k + r 



(-l)'27r(A; + l) 



_ 7r'^{k+iy + {'&-^'y ^7r2(A; + 2A)2 + (^?-?9')2 

_ — k — 



^ OO 

^E 

k=0 



k+1 



+ 



1 



-1^ 



tf-^' Y k + 2A 

K(k + l) J 



■d--d' 

7r(fc+2A) 



(4.24) 



Some comments are to be made. First of all, even powers of j do not appear in this 
expansion, because the energy density is written as the sum of two odd function. Of 
course, it follows that all even powers in the interacting case will vanish when we shut 
down the interaction. What is more, the first three term are somehow "stable" (see 
comments to fl5.42p against the insertion of a less trivial S-matrix. 



5 The very first orders in B 

We recall that in the 0(6) strong-coupling regime we have B «1. This way written, 
the kernel allows a easier calculation of the first few orders in B of the quantity s{'d). 
The first, raw approximation consists in writing the solution e{'d) as composed of the 
forcing therm alone, e{'d) = h — mcosh'i? The boundary condition implies 

B = arccosh ( -] ^ J 2^?—^ (5.25) 
\m J \ m 

This bare result is already enough to obtain the fermion gas approximation, as already 
stated. We would like to know how much this result for the energy density is reliable, 
i.e. if calculations involving higher orders in B are likely to leave the first coefficient as 
they are, or else if they are destined to upset the solution. We will go on by brute-force 
Taylor expansion in B. To begin, we illustrate the second order case, an then we will 
give the systematics of the method. 
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5.1 What we can do up to the order 

By defining the new variable x = ^ and expanding the fractions that compose the 
kernel in B we obtain, up to the order 



/i — m cosh + i? J K{d — Bx)e{Bx)dx = 

h — m cosh + B J K{{} — Bx) — mcoshBx + B J K{Bx — By)e{By)dy^ dx 



h — m cosh H — - / dxy^ 

^ k=0 

( ^ „i oo 

h — m cosh Bx H — - 

B °° 

h — mcos\l'^^ — ^> / 

u—n — i 



k+ 1 



. I ( ^-Bx Y k + 2A f^^B^Y 



1 



k+ 1 



, f Bx-By \ 
^ ^ I, 7r(fc+l) ) 



k + 2A 



, f Bx~By_\ 
^ ^ \^7r(fc+2A) ) 



\ 



e{By)dy 



dx 



J 



-If 



k+ 1 



d- Bx 
Tv{k + l) 



+ 



d-Bx 
n{k + 2A) 



'If 



k + 2A 



e{By)dy ) dx 



h — m cosh i3x H r- > / 

/i - mcosh^? + - / \Si-SY'&- Bxf] 

B 

h — m cosh(_Bx) H / Si{h — m cosh Sy) dy 



(5.: 



where 



Si = -5: 

k=0 



(-ir , (-ir 



A; + 1 A; + 2A 



and 



£(■(9) ^ /i — mcosh-i^H- 
5 .1 



h-m{l + 



fc=0 



+ 



Si {h — m) dy 



2 



j ^ [Si -SY'd- Bxf] 
B 

h - mcosh^9 + - / ISi-Ssi^- Bxf] 
h - mcosh^? + - I 2Si (/i - m) - — ^i^^ /" x^rfx + 4-3? (h-m) + 

TT y 2 J_i TT 



(/i - m) + 2— Si {h - m) 

2 TV 



dx 



~S3{h-m) I {^-Bxfdx + ^S^B^ J {d - Bxf x^dx + 



B \ 

-2-S1S3J {^-Bxfjdx 



(5.27) 



1 „o\ 2 . „ „ 4 



e{B) = h-mil + ^B^] +-{h-m)SiB + —{h-m)SfB^ + o{B''^) 

\ 2 J TV TT^ 

= h-m + -{h-m)SiB+i-— + —{h-m)SnB^ + o(B^) (5.28) 

71 \ 2 TT^ J 

By solving the condition e{B) = we gain an expression for the extreme 



-{h -m)^ + sJ^Sl{h - my + f (/i - m) - ^{h - mySl 

2 



3 ^ 



2., X 2^/i-m ^ fh-m\^. 

-ih-m) + -Si + 0\{ (5.29) 



m n m \ \ m 



where we have coherently omitted all contributions from B^ on, that is (^^^p) ^ and 
higher powers. For what the free energy is concerned, we can write the BA formula and 
approximate to the second order in B 

B 1 

f (h,m) = / cosh{}e (■&) d-& = B cosh (B xq) e (B xq) dxo 

m J J 

} ( f \ 

B cosh(i?a;o) I /i — m cosh (5 Xq) + -B / Si{h — m)dxi dxo 

~ B + i^-m- ^B'^xo^ + B j Si{h - m)dx^ dxo 

1 

^ Bj{h-m + 2BSi{h - m)) dxo = 2{h - m)B + A{h - m)Si(B'^0) 
-1 

the density of thermodynamics, up to the same order, is 

B 1 

j = —^f{h,m) = ^ J d^e {^) coshi? = ^B J \ 1 + SiB j dy \ dx 

-B 1 



{B + 2SiB'^) (5.31) 



m 

SO that the energy density becomes 
1 / 1 



n = J ll + SiB J dy\ -^{2{h-m)B + 4{h-m)SiB^) +0{B'^) 



m 

7T 



1 

2 



{B + 2SiB^) + O {B'^) = mj + 0{f) (5.32) 
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Indeed, contibutions of order are all included in the mj term. Unfortunately, to 
catch the contribution we have to keep the B^s, as well as to catch the contribution 
we have to keep the ^"'s, an consider at least (n-1) nested kernels in the pseudoenergy. 

5.2 What we can do up to the order 

Once again, we briefly repeat the very same steps as before, reaching one higher power 
of 5. 



£ = /i— m cosh -(9+1? J dxK — Bx) 

-1 



h — m cosh [Bx] + B j dyK [Bx — By) [h — m) 

(5.33) 



/ ^ ^HLb [dx(l- -B'x' 



h-m- ^B'^x^ + B j dxK {Bx - By) 



771 1 

h-m B^x"^ + 2B-Si {h - m) 



TT 



1 

^ ^HLb [dx(l- -B'x' 
2n I V 2 



771 1 /I 

h-m B'^x^ + 2B-Si {h-m)+ AB^ -Si ) (h - ng> 



TT 



TT 



34) 



Of TfL f 

j = —-T— ^ — / rf-i^cosht? 
dh 271 J 



1 + 2B-Si + AB^ ( -Si 

71 V TT 



B + -B'^ 
6 



1 + 2B-Si + ( -Si 

TT \ TT 



m 

7T 

-(b + 2-SiB^ + A^S^B^ + -B' 
n \ 71 71^ 6 



(5.35) 



1 

n = f + jh- —B dxil- -B^x^ 

-1 

2 

m o 
= mj + —B^ 

OTT 



1 + -B^x^ + 2B-S1 + AB^ i-Si 

2 71 \7l 



(5.36) 



now, the term with B^ can come only from the third power of the density j, because 
we have already acknowledged the term linear in j, so we calculate the coefficient of 

2 

the j-expansion of as the ratio among ^ and the coefficient of the first power of B 
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in the expansion of ?': fiq = ^ 



So we have 



TT 



= mj + — j3 + O(j^) 



(5.37) 



5.3 A small improvement 



We saw in ( l5:29|) that B ^ + O (^). We can use this to show the stabihty of 

the non-interacting fermion gas approximation and that the most crude approximation 
yelds a correct coefficient of j^. 



e{d) = h-m cosh + j d^'K - d') h-m cosh +B j d§"K {§' - §") {h - m) 

-B . -1 



h — m cosh d + B I d-d'— Si 

TT 



h-m cosh {'d') + 2B-Si {h - m) 

TT 



+ 0(B^ 



h — mcosh^} + B—Si / dx 

TX 

h — m cosh d + R{B) 



h — m cosh [Bx) + -B^i / dy {h — m 



0{B^ 



(5.38) 



where R{B) = 2 {hB -msinhB) + AB^ {h-m) ~ 2(/i - m)(i? + 2B'^) behaves as a 
modification of the chemical potential. Thus, we can expect that our result will not be 
different from the free one. Indeed 



B 



fih) 



Tfl f 

/ d'd {h + R (B) — m cosh {}) cosh {} 

27r J 



B 



B 



m 



27r 



— / d^\h + 2{h- m){B + 2B') -ml B' 1 B' + 0{B 



^2 



^2 



777 

- {{h - m)B + 2{h - m){B + 2B^)) + 0{B^) 



We easily calculate, from (14.221) . 



(^9) = 1 + 2SiB + ASiB^ + 0{B^) 

oh 



and, from fl212ll 



(5.39) 



(5.40) 



B 



(J 771 C Cf 771 

--^f = — U^e (d) coshddd = -Bil + 2SiB + ASiB^) + O(B^) (5.41) 
oh 2% J oh TT ^ 
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to conclude from (12.131) . somehow more "efficiently", that 

777 

n = jm {{h - m)B + 2{h - m){B + 2B^) + {h - m) 2B (l + 2SiB + 4:SiB^))+0{B^) 

(5.42) 

now, by remembering that h — mc:^ mB^ + 0{B^), it becomes evident how the structure 
of the forcing therm ensures the stability of the free-fermion approximation: every 
correction affects the series from the -B^ (i.e. from the j^) term on. Moreover, due to 
the absence of the Bs and of the -B^s, it suffices the ffist-order approximation for B to 
fix the coefficient of j^. 



6 A method for calculating B to all orders 

We have seen that we can express B in powers of x = ^^^p as 

oo 

5 = ^6„x" (6.43) 

n=0 

as well as we can write the pseudoenergy as a power series 

oo 

e(5)=m^e„5" (6.44) 

n=0 

At the same time, the x-dependence of the coefficients can be explicited by Taylor 
expansion. This means that 

e {B{x)) = e[,°) + el^^x + e'^^x^ + ...+ (ef^ + e^^^x + ef^x^ + . . .) (6o + hx + hx^ + ...) + 
+ (ef + e'i^x + efx^ + . . .) (6o + Kx + h^x^ + ...)' + ... (6.45) 
We saw that, up to the second order, the only nonzero coefficients were 

.(2)_i p(2)_2c p(0)__l .(2)_Aq2 

SO we can solve, order by order, the boundary condition (12.81) 

zero ~i^o = ^ 6o = 

one = useless 

two e^o^ + e?bi + efb2 + ePbl + bo{...) = l-lbl = ^ bi = V2 

three ef^bi + 2ef'>bib2 + efhl + bo{...) b2 = 2^ + 2ef 

the zeroth order fixes bo: this guarantees that it is possible to calculate bk only with 
the ffist k orders in B of (12.81) . In general 



1 " 

&n-i = --^5^E(l-^™Ao)e(^) J2 (6-46) 



2e^\ 



2 "lm=li=0,2 jj+...+j„=n-t 
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with 

p(0) 
.(2) 



7r"(2»)! 

p=0 n+2p+2A:i + ...+2fc„=Ar ^ ^' 



dXi... I dXn (1 - Xif^^ . . . {Xn-1 - Xn)^"" X^^ 



\2fc„ t'^P]^^ 



E 

n+2fci+...+2fc,i=7V 



TT'' 



dxi . . . y dx„ (1 - xO''^ . . . (a:„_i - x„)^'^" fi^^ (6.47) 
1 -1 



7 ^(j)' a systematics to all orders 

By substituting the exphcit formula for the kernel in (I4.19P 



^ CXD 



k=0 



A; + 1 



^ ^ l^7r(fc+l) y ^ l^7r(fc+2A) J 



7r(A; + 2A) 



(7.48) 



we get 

e{'d) = h — mcosh.'d + 



i 1 

OO -| « oo 

+ dx,... dx^Yl 

n=l 1 "^-(^ mi=0 



E 

m„=l 



1 



m„ + 1 



1 



7r(fe+l) 
1 



A; + 1 



+ 



1 



nir, + 2A 



1 + 



[h — m cosh {Bxn j) 



/i-mcosh^9 + ^— 5" / dxi... / ^ (-1)"^ ^ (-1)^'^ 



n=l 
1 / ?9 - 



mi + 1 \7r(mi + 1) 



d-d' 



mi + 1 \7r(mi + 1) 
and by defyining the series 

oo 

s» = E 



1 

2ii 



2ji 



_1 -mi=0 ji=0 

1 / d-d' V'' 

mi + 2A V7r(mi + 2A) / 



d-d' 



mi + 2A V7r(mi + 2A) 



2ii 



. mi=0 ji=0 

(/i - m cosh (5a;„)|7.49) 



fc=0 



+ 



TT*^ (A; + 1)" TT" (A; + 2A)' 
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we can rewrite the previous result in a much more tangled way, but with a more 
transparent insight on the contributions of the different powers of B 

e{'d) = h — mcosh.'d + 

n=l ^ ji-jn=0 ^ ^ 

. . . S,„B'" {xn-i - Xnf'" {h-ni cosh {BXn)) (7. 

Please note that we are allowed to exchange the order integrations and summations 
over the ji indexes, because on the finite interval [—1, 1] all powers of x are integrable 
and the succession of reduced sum does converge to the actual kernel. 

s{'d) = h — mcosh.i) + 

oo oo } } / l') \'^^^ 

+ E- E i-iy"B^+''-S2,,^,S2,,+i...S2,^+i dx,... dxJ--BxA ... 

n=l ^ h-iu=0 11^ ^ 



2p 



...{Xn-l-Xnf^-ih-mY^B'^.^ 



h — m cosh d + 



(-1) " f f fd 

+ E E ^n-j„ S2j^+iS2n+i--S2j^+i I dxi... / dxn {X1-X2) 

n=lh...j„=0 ^ " _i _i ^ ^ 



. . . - Xnf- \h - m j (7.51^ 

where Jn = Ji + J2 + ••• + jn- In the following we will also define 

Sn{{j}) = 5'2ji+l5'2j2+l---5'2j„+l (7.52) 
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following [19] we can compute the free energy as a function of the chemical potential h 
and of the mass m of the field 



2tx 
m 



f (/i, m) = J cosh -i^e (-i?) di!} = B J cosh (B xq) s {B xq) dxo = 

-B -1 

{oo ^ oo oo \ \ 

^Et^E^ E i-iV'^sAij}) dxo... dx^ 
r=0 ^ n=0 jij2-jn =0 

...X^J (Xo - X,f^' . . . - X^f'^ 5l+n+2J„+2. ^ 

co^oo-^oo ^ \ \ \ 

-^Et^E^ E i-^y"Snm)Y.7^J d^o... dx^ 

now it is to the density j, dual to the chemical potential h through the Helmholtz free 
energy /. 

B B / B 



j dde {d) coshd = ^ j dd{l + j dd'K - d') j cosh^9 



-B -B \ -B 

11 



SEt^E;^ E i-iy^sAij}) dxo... dx. 

r=0 ^ >' n=0 iii2...in =0 



We have used ( I4.23P and (12.81) . The Legendre transform follows straightforwardly 
^{hrn) = f {h,m) + jh = 

1 1 

2 00-.oo-.oo oo-.« « 

r=0 ^ ^ n=0 jii2...in =0 p=0 ^ ^' 

Xl' (Xo - Xi)^^-^ ■ ■ ■ - Xnf'^ 4P5l4-n+2.„+2.+2p ^ 

m^°°l°°l°° ^l} } 

r=o ^ ^ n=o iii2-in =0 p=i ^ -f^^ 

X'J (Xo - Xi)'^'^ ■ ■ ■ {Xn-1 - XO'^'^ ^2p^l+n+2J„+2.+2p ^7 55) 

Please note that the p index in the last summation, having recognized the p = term 
as the linear contribution in j, now runs from 1 to infinity. This, of course, allows us to 
exclude any term of order j'^, as the lowest power of B in the summation is now three. 
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At this step, we have the two series 

oo 

n.=l 
oo 



(7.56) 
(7.57) 



n=l 



In this situation, one would naturally appeal to Lagrange inversion formula and express 
Q in powers of j as 



d 



n—l 



(7.58) 



B=0 



n=l n=l 

where 0(5) = 

7.1 Some explicit coefficient 

However, in order to perform some calculation, it may be simpler to extract the coeffi- 
cient of the n-th power of j, from the coefficient of the powers of B up to the n-th: 
uji, . . . , ujn- In fact, having already isolated the term mj, we easily calculate fl^: 



uJi = ^ I dx 

Ztt , 
-1 







m 



TX 



^3 



vnr 1 
2^2 



xo'^dxQ 



m 



We can make another step without too much trouble 

2iT 2 IT / 67r^ Sn'^ 



(7.59) 

(7.60) 
(7.61) 



m 



Si 



(7.62) 



^ n IT 2% IT J J 71^ 

-1 -1 



(7.63) 



we will now specialize the solution to the 0(6) sigma model, i.e. we set = 6, so that 



A = i 



^ oo 

vr 

fc=0 



k+1 



1 , 1 
- n2 + - 
vr 2 



(7.64) 
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^4 - ^3(j^)4 1 



Jl 



m 
3^ 



.Si 



6m vr^ 



5i 



2^ 
3 m^ 



1 , 1 
- n2 + - 

vr 2 



(7.65) 



where we have written the coefficient of in the series in B of as {j^)4- Now we 
have warmed up, we can approach the next one: 



CJ5 

if). 



— (— + 
TT V 120 3 7r2 ^ 



2 1 

3^ 
m\3 m 



TX 



TX TT 



-Si 



/m\3 /12 



12 

TT 



— C2 



(7.66) 



— (cj5 - ^n-ijlja - sn^jij^ - 4^4^52) = 

Jl 

(7.67) 

On one hand, please note that the coefficient Q5, as well as those of all odd powers of 
j, does reduce to its free approximation if we switch off the interaction, i.e., if we send 
S*! — > 0. On the other hand, since the oddness of the non- interacting gas series, the 
Q4, correcly vanishes in this limit. Generally speaking, unfortunately, when considering 
the order B^, we need to subtract all the terms coming from lower powers of j, which 
makes the calculation quite cumbersome. 



^5 = -5(^5-fi4(/),-^^3(/)5) = 

m3 I 40 7r2 V 



7.2 Comparison with Lagrange formula 

Of course, this results does match with those coming from Lagrange formula. The ffist 
coefficient is trivial, the second reads 



fig 



1 d 
2dB 



P=i 

' 00 



1 
2 



J 5=0 



p=2 



,P=1 



k=2 



vfc=l 



vfc=l 



1 / 2ct;2 ?2 , 
' ^ -2—] = Q 



2 V Jl Jl 
Checks. The third 

1^3 = 



(7.69) 
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O'(0) = uoi 
O"(0) = 2cu2 



n\UJn 



(7.70) 



$'(0) 
$"(0) 



d B 

■ 

Jl 



(i)' 



_B=0 



J2_ 



(7.71) 



6— 3 

Jl 



12 



Jl 



3^1 ( 



,J2 -Jlj3 
' Jl^ 



J B=0 



^^3 - J3 



,■3 
Jl 



12 



Jl 



12 



1 TT^ 



6 rrc 



m 2 



U^J2_ 
■4 

Jl jr 

m 1 



27r3 271 7r 



271 3 271 7l2 



checks again. We will not go further. 
7.3 Analiticity 

The convergence of the two series in powers of B will be estimated by combinatorial 
considerations. We have possibilities to arrange k nonnegative integers in such a way 
that their sum yelds N . Moreover, the number n of the factors appears in the exponent 
together with the factors j„ themselves and the two indices p, s. So the general A^-th 
coefficient grows less than 



N 

o2n 

n.=0 



2^" {N - n) 



n+2 



where we have used 



1 


1 


j dxo . 


■/ 


-1 


-1 



and 



Sn 



dXn-- ■ {Xq - Xi 



E 

fc=0 



>2ii 



(a:„_i - xj'" < 2" ■ 2 



71" (A; +1)" 7i"(fc + 2A)" 



< 1 
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We find the maximum value of tlie addend 



— (N-x 
ax 



,x+2 



fix" 



{N - xf 



x+2 



e(^+2)i°(^---)(ln(iV-x) 
d 



x + 2 



dx 

(x+2)\n{N-x) 



{x+2)MN-x)l^^^^_^-^_ 



N -x^ 
x + 2 



[N — a) = e^-a 



N-x 



\n{N - x) 

g{a+2)ln(JV-a) (iV - a 



X + 2 



N-x 
2 



2 ^ x + 2 



N-x {N-xY 



+ -T^^^ ) < (in X = 47.73) 



N-a {N-xY 



and we substitute it in tlie sum 

N N 



(AT 



n=0 



n=0 



^ ^2N ^ < 22^ (iV + 1) e^-;^ ~ 



N{Ae 



(7.74) 

Tfie general term grows witli tlie N-th power, so that the two series of powers of i? ^ 1 
easily converge. To conclude, Lagrange theorem expresses the analiticity of the series 
f l7.58p . when written in powers of — <C 1. 



8 Checks and previsions for strong SYM4 

As noticed by [1], the strong coupling limit j << g and the consequent reduction to 
the 0(6) bosonic nonlinear sigma model in two dimensions allows to calculate exactly 
the the anomalous dimension of high spin operators f{g,j) with j = j^'- the S0(6) 
charge density modifies the scaling by a 2Q, 

fig,j) = Mg) + 2n{g,j) (8.75) 

The possibility of calculating, at least in principle, many successive coefficient /i, /2, • • • 
of the expansion (as announced in the footnote 10 of [2]) allows a check of the first 
scaling functions from the SYM side, appearing in the Sudakov scaling (11. 4p . At this 
purpose, we underline that, on the one hand, the density j that we have used is defined 
from the string side in (3.1) of [T] as 

^ (8.76) 



2 log 5 



while on the other hand, recently available calculations [3] approach the problem with 
slightly different notations. From their point of view 

jsYM = = 2j (8.77) 
log^ 
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It follows that, rescaling their results, we get an explicit link with the 0(6) model, 
namely 

/„ = 2^-'n^ (8.78) 

The first three coefficients from the sigma model point of view (fii, Q2 and ^3), are 
exactly the same as for the free, non-relativistic theory. They were already given as an 
approximation in [19], as we said above. Further, we have performed the calculation 
of Q4 and Q^, that contains trace of our interaction and the model-dependence, too, 
entering through 5*1 as in (17.641) . On the other front, the correspondence with /i, /2, 
/a and /4 is a remarkable fact. In particular, this very last scaling function catches the 
nature of the system and selects the specific model. 
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